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Development of equations for predicting fuel bed depth (called ''bulk 
depth" herein) appropriate for modeling fire behavior in slash is de- 
scribed. Bulk depth (y) was correlated with the expected number of 1/4- 
to 1-inch-diameter particle intercepts per foot of vertical plane transect 
(x) by regressions of the form y = avx. Values of "a" suitable for use in 
fire models were 0.767 for high-lead harvest debris, 0.940 for precom- 
mercial thinning of pines, 1.22 for precommercial thinning in several 
other western conifers, 0.877 for ground-lead harvest debris in pines, 
and 0.542 for ground-lead harvest in other species. Lopping of slash 
reduced average depth 17 percent for harvest debris and 31 percent for 
precommercial thinning debris. Correlation of high intercept depth 
(maximum height of sampled fuel particles)with bulk depth showed that 
the bulk depth can be well predicted using 64 percent of the more easily 
measured high intercept depth. 

Models for settling of slash, retention of foliage and fine twigs, and 
species mixing, useful in preparing data for fire models, are presented. 
Application of the models in a slash hazard appraisal computer program 
is illustrated. 


INTRODUCTION 


Slash or debris created by harvesting and thinning are a major fire management 
problem because these residues can create unacceptable fire behavior hazards. Treating 
slash to maintain an acceptable fire hazard is expensive and requires skillful decision- 
making. An inexpensive, simple-to-use, yet objective means of appraising the potential 
fire behavior of slash is needed to aid decisions in managing slash. Knowledge of 
potential fire behavior can help determine treatment alternatives, the financing of 
slash treatment activities, and even determine whether the slash should be created. 

This report describes a method for predicting depth of slash fuels for analytical model- 
ing of fire behavior. 


The capability to predict debris and to model fire behavior has made possible a 
quantitative system for appraising fire behavior potential. Rate of fire spread and 
flame front intensity can be mathematically modeled using Rothermel's (1972) fire model, 
which is the basis for computing spread and energy release indexes in the National Fire- 
Danger Rating System (NFDRS). An area-growth-rate model by Anderson,! crown scorch 
model by Van Wagner (1973), and flame length models by Byram (1959) and Thomas (1962) 
are also available. A method for assessing total heat release and period of flaming 
that incorporates large-diameter fuels has been theorized by Albini.* The fire models 
require fuel loading by size class and fuel bed depth as inputs. 


Loading and depth of debris vary considerably by cutting prescriptions, therefore 
must be estimated for each individual cutting situation. Other fuel inputs such as 
particle density, heat content, silica-free ash content, and particle surface area-to- 
volume ratio tend to be species dependent and can be approximated from known values 
(Albini 1976; Brown 1970b, 1974). 


lOffice memo of August 10, 1973, on file at the Northern Forest Fire Laboratory, 
Missoula, Montana. (Manuscript in preparation by Hal E. Anderson.) 

2Albini (1976) outlines the total heat load computation; the computation of burning 
time is based on empirical work reported by Harmathy (1972) describing structure fires 
and pile burning, modified slightly to agree with single-particle burning times reported 
by Anderson (1969). 


For Rocky Mountain conifers, loadings of slash can be predicted from relationships 
between tree crown weight and tree characteristics such as d.b.h. and height (Brown 1978). 
In the USDA Forest Service Northern Region these relationships have been installed in a 
computer program that obtains input from tree inventories and predicts potential debris 
as output.°’* In addition to predicting debris, a method to predict fuel depth is needed 
to appraise fire potential. The objectives of this study were: 


1. Determine the relationship between fuel depth and loading of slash and the 
extent to which species, age of slash, method of skidding, lopping, and other factors 
influence the relationships. 


2. Determine the relationships between the easily measured high intercept depth 
and bulk depth. 


When using analytical models such as Rothermel's (1972), fuel depth is a critical 
parameter because it determines bulk density of the fuel array for given fuel loadings. 
Rate of fire spread is very sensitive to bulk density (Williams 1977). Fuel depth is a 
measure of the vertical extent of fuel in the zone that is actively involved in the 
spreading flame front. Conceptually, the bottom of this zone is the forest floor and 
the top is the height where fuel ceases to exist or is too sparse to affect propagation 
of the flame front. Fuel depth that is compatible with fire modeling can be difficult 
to measure because locating the top of this hypothetical zone in the fuel array requires 
judgment. 


Fuel depth in slash and other downed woody material has been measured primarily on 
a high intercept basis (Brown 1974). In this procedure, the top of the fuel is defined 
by the highest particle to intersect a vertical plane about 1 foot wide. Although this 
procedure is easy to learn and to use, it does include large void spaces in the fuel 
array whenever they occur. This permits overestimation of an effective fuel depth and 
underestimation of bulk density required for fire modeling. A procedure for measuring 
effective fuel depth--called "bulk depth''--has been developed by William Frandsen (1974) 
at the Northern Forest Fire Laboratory. It allows observers to account for void spaces 
in the fuel arrays and provides a measure of effective fuel depth appropriate for fire 
modeling. 


The merchantable top diameter relates to depth. Small top diameters result in more 
lopping of supporting branches and removal of bolewood than large top diameters; hence, 
the slash is more compacted. Species may influence depth of slash due to differences 
in branch stiffness and branching habit. In a study by Roussopoulos and Johnson (1975), 
loading and depth of slash were directly related; however, as loading increased, depth 
increased at a reduced rate. Probably as more tree crowns are added to a fixed area, 
overlapping of branches occurs and the weight may cause compression. Methods of felling 
and skidding trees should influence depth because the amount of trampling and breakage 
depends on if and how merchantable boles are removed. 


Settling with age has a most significant influence on slash depth. For some 
western conifers over a 5-year period, depth of lopped and unlopped slash was reduced 
to one-half of the original value (Fahnestock and Dieterich 1962; Kiil 1968). Rate of 
settling varied; for some species, depth actually increased slightly during the second 
year before settling. In a study of piled slash over a period of 29 years, Wagener 
and Offord (1972) found that piles continued to settle; however, 50 percent of the 
settling occurred during the first 5 years. 


3Users' guide to debris prediction and slash hazard appraisal. 1977. USDA Forest 
Service Northern Region, Division of Fire Management, Missoula, Mont. 

*Brown, James K., and Cameron M. Johnston. 1976. Debris prediction system. Office 
report on file at the Northern Forest Fire Laboratory, Drawer G, Missoula, Mont. 


METHODS 


Fieldwork 


Sampling destgn.--The study was designed to provide data for four different skid- 
ding methods: 


Ground lead.--Crawler tractor and rubber-tired skidder, entire log dragged. 
Skyline.--One end of log elevated. 

Heltcopter.--Limited skidding. 

Precommeretal thinntng.--No skidding. 


BWNH 


For each skidding method, two species groups were desired and for each species 
group, two age classes of slash. For each combination of skidding method, species, and 
age, we attempted to locate two study sites on areas represented by each of three aver- 
age loading levels: low (sparsely distributed), medium (nominally about half of ground 
covered by slash), and heavy (uniformly distributed). 


Study sites were located in areas having different loadings to assure that any 
differences in skidding patterns due to loading would be reflected in the data. 


Species were partitioned into the following two groups, based on similarity of 
lopped depths as observed by Fahnestock and Dieterich (1962): 


1. Douglas-fir (Pseudotsuga menztesti (Mirb.) Franco) 
Grand fir (Abtes grandts (Dougl.) Lindl.) 
Alpine fir (Abtes lastocarpa (Hook.) Nutt.) 
Western redcedar (Thuja pltcata Donn) 
Western hemlock (Tsuga heterophylla (Raf.) Sarg.) 
Western larch (Lartx oeetdentalts Nutt.) 
Engelmann spruce (Picea engelmannti Parry) 


iS) 


Lodgepole pine (Pinus contorta Doug.) 
Ponderosa pine (Pinus ponderosa Laws.) 
Western white pine (Pinus monttcola Dougl.) 


All slash was at least one-half year old and had existed through at least part of 
a winter. Two age classes were recognized: O- to l-year and 3- to 4-year. Some combin- 
ations of skidding, species, and age were not found. Because skyline and helicopter 
logging has received substantial use only recently in the Intermountain area, the O0- to 
l-year and 3- to 4-year lodgepole pine and 3- to 4-year Douglas-fir were unavailable for 
sampling. Surprisingly, 0- to l-year precommercially thinned stands other than pines 
were also unavailable. All other skidding, species, and age combinations were sampled. 


Within study areas that were a minimum‘of several acres in size, two primary tran- 
sects were established. Along each primary transect, a total of 50 or more sampling 
points were located at 2-foot intervals. Points without slash were not sampled. 


Loading measurements.--Using the planar intersect technique (Brown 1974; Brown and 
Roussopoulos 1974), loading was measured for at least 15 randomly preselected sample 
points along each primary transect. Sometimes the random selection resulted in more 
than 15 points out of 50 for measurement of loading. For 1/4- to l-inch particles, two 
2-foot planes, crossed perpendicularly, were vertically oriented and intersections 
counted. For l- to 3-inch and greater than 3-inch pieces, the 2-foot planes were ex- 
tended to 4 feet. Particles 0 to 1/4 inch were not tallied because of the counting work 
involved and the fact that the information appeared unnecessary to meet objectives. 


Depth measurements.--High intercept and bulk depth measurements were recorded at 
each sampling point. High intercept depth was measured as the vertical distance from 
the bottom of the litter layer to the highest O- to 3-inch diameter slash particle 
intersecting each 2-foot plane. Pieces greater than 3 inches in diameter were omitted 
in determining depth because they occur infrequently as the highest particle and have 
considerably less influence on rate of spread than smaller pieces. 


Bulk depth was measured in each of four pie-shaped quadrants of a 2-foot diameter 
cylinder whose central axis was vertically oriented at each sample point. The two 
perpendicular sampling planes for tallying intersections of 1- to 3-inch particles, 
delineated the cylinder into quadrants. The top of fuel was the average height of an 
imaginary pliable sheet draped over the fuel particles. The bottom was at the base of 
the litter layer. Vertical gaps free of fuel for more than 1 foot were subtracted from 
each quadrant's depth. Gaps of less than 1 foot were assumed to maintain vertical con- 
tinuity of flames, thus were included in the depth measurements. Depths of the four 
quadrants were averaged to obtain a bulk depth estimate for each sample point. 


Lopping.--After loading and depth were measured initially, the slash along the 
transect was lopped so that all branches were within 2 feet of the ground and boles 
within 1 foot. High intercept and bulk depths were remeasured at all sample points 
where depth had changed. 


Analysis 
Rattonale 


To understand the rationale behind the method of analysis applied here, bear in 
mind the objectives of the effort and the nature of the variables. The first objective 
was to establish an equation that can be used to predict the mean bulk depth of a slash 
fuel bed from quantities that describe the amount of slash on the area. The sampling 
procedure provided an estimate of the bulk depth and the variables for quantifying the 
slash loading at each sample point: the number of intercepts of 1/4- to 1l-inch and l1- 
to 3-inch-diameter fuel particles. The expected number of such intercepts can be 
predicted from slash loading, tree species, and d.b.h. simply by parsing the loading 
according to the fractional weight distribution of the individual tree crowns (Brown 
1978). But the number of intercepts (regardless of size class) is subject to great 
sampling variability, since it should be approximately Poisson distributed. 


A Poisson-distributed variable has a variance equal to its mean, so if N is its 
expected value, samples from N - YN through N + YN should occur with approximately equal 
frequency. This intrinsic variability makes it impossible to distinguish with certainty 
between a change in the mean value of the number of intercepts and simple data scatter 
when moving from one sample point to the next. To alleviate this problem, the data can 
be aggregated, combining measurements that lie within a range +VN of each other. We can 
use the average values to discern the underlying trend of the average bulk depth with 
the average number of intercepts; otherwise, the trend would be largely obscured due to 
the great scatter in the individual samples. Sampled bulk depths also exhibited 
Significant scatter, indicating a need to aggregate data. 


Fuel bed depth was related to the expected number of intercepts of 1/4- to 1l-inch- 
diameter fuel particles per foot of planar intersect by fitting regressions through 
aggregated data points. The 1/4- to l-inch particles served as a proxy for the loading 
of debris fuel under 3 inches in diameter. The-fitting of data points resulted from a 
three-step averaging process as described below: 


1. For each study area, the data from the two transects were treated as a single 
set. Only sample points at which intercepts were counted were assembled in the data 
set, so the "unit" of sample information consisted of a triplet of numbers--the bulk 
depth, in inches; the total numbers of intercepts of 1/4- to l-inch size class fuels in 
the two crossed sample planes; and the number of intercepts of 1- to 3-inch size class 
fuels in the same two (extended) planes. 


For inspecting these triplets of numbers, they were aggregated and displayed in 
five categories, according to the number of 1/4- to l-inch intercepts: 0-2, 3-6, 7-13, 
4-22, and 23 or more intercepts. The display (table 1) consisted of the following 
descriptors for each subset: 


Average bulk depth 

Average 1/4- to l-inch intercept count 
Number of sample triplets in the subset 
Mean square bulk depth 

Mean square 1/4- to 1l-inch intercept count. 
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Scrutiny of such tables quickly revealed sample points that were obviously not rep- 
resentative, so they could be culled from the data set. The mean square depth figures 
served well to permit the ''outlaw'' points to be identified. Fewer than 10 points were 
discarded from more than 1,500 collected. 


2. These data were analyzed in many ways in attempting to discover trends and 
correlations. One fact that soon became evident was that no substantial dependence of 
bulk depth on the 1- to 3-intercept count could be established, whether or not the 1/4- 
to l-intercept count was included. This allowed further simplification by combining 
data for all 1- to 3-intercept counts. Table 2 displays this simplified data aggrega- 
tion. An asterisk indicates that a sample point has been discarded in computing the 
averages shown and two asterisks indicate the discarding of two data points. These data 
form the basis for the regression relationships between bulk depth and number of 1/4- to 
l-inch intercepts. 


3. Data in table 2 were grouped by combinations of skidding method, species, and 
age of slash. Constrained regressions were applied to the data because logically if no 
1/4- to l-inch fuel particles exist there should be no bulk depth. Analysis using un- 
constrained regressions supported this approach in that the constant terms were small 
and often statistically nonsignificant. 


In deriving the regression equations, the data from all units aggregated were 
averaged. That is, all the bulk depth-intercept count pairs in the intercept count 
range 0-2 were averaged, all those in the range 3-6 were averaged, etc. Then the re- 
gressions were formed by weighting each resultant aggregate data point by the number 
of its supporting measurements. 


Scatter diagrams of these "grand average'' points showed a fraction power law form, 
so we chose to regress the average bulk depth against the square root of the average 
number of 1/4- to l-inch intercepts. 


1/4 to 1 
intercept 
count 
range 


Table 1.--Sample table of aggregated data from one area used 
for establtshtng the relationships between bulk depth 
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TABLE 2.--Results of partitioning depth and intercept count data by intercept count range. WN = nwnber of data points, X = average count of 1/4- to 1-inch 
intercepts, D = average bulk depth (inches), D* = mean square bulk depth (inches*). Scanple wits are grouped on the basis of similar types of 
cutting, species, and debris age. An asterisk, *, denotes deletion of one measurement point, ** denotes two points omitted 
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1975-11 
1975-13 
1975-14 


1975- 7 
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1975-12 
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11 4.36 5.58 57.66 14 9.36 4.32 23.48 4 15,25 6.78 48.43 
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11 4,18 2.70 9515, 9 8.11 5.92 40.90 5 15.60 12,88 173.7 
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6 4.17 Se07: 14.74 25 8.76 19.32 439.7 2 14.50 21.25 469.6 

5 4.80 12.42 208.9 18 9.83 14.62 266.6 13 16.69 12.02 186.4 2 23.00 14.15 202.0 
8 4.50 6.85 75.78 12 8.83 13,12 204.2 1 14,00 13,80 190.4 1 32,00 S.80 33.64 
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a S25: 7.90 90592 11 10,45 18.35 430.6 8 18,13 22,10 552.3 4 28.75 24.10 614.7 
13° 4.38 6.05 58.07 7 9.57 11,57 188.4 S 19,00 13.08 185.1 4 28.25 21.18 457.9 
8 5.25 7.13 85.75 3 9,67 19,00 369.7 


One penalty for using constrained regressions is that there is no widely accepted 
simple quantity that reflects the degree of agreement between the equation and the data, 
in the way that the "coefficient of determination," called r*, does for unconstrained 
regressions. If Ye represents a bulk depth data point to be fitted by regression, Wears 
the value of the bulk depth predicted by the regression expression, and y the average 
of the data points to be fitted, then, by statistical theory the expression for the 
coefficient of determination for a constrained regression should be 


Die 7 ee. 2 


This quantity is, however, a poor measure of the "goodness" of the regression equation. 
So we have used, in this presentation, a nonrigorous but intuitively more appealing 
measure of the suitability of the regression description: 


S.= us Bs Va G/ Or ae le 


This parameter compares the variance about the regression to the variance about the mean, 
as does "r2" for an unconstrained regression, but it is not limited numerically to the 
range 0 to 1. It is restricted only to be less than unity, and can be negative. 


RESULTS AND DISCUSSION 


Depth Versus Intercept Counts 


Regression results for all of the skidding, species, and age combinations obtained 
in the study are shown in table 3. The relationships for skyline and helicopter 
skidding were similar; so they were combined and called "high-lead.'"' The high values 
of the fit descriptor, s, are largely due to the manner of aggregating data. Examples 
of fit are shown in figures 1, 2, and 3. Because of the data aggregation, tests for 
differences among the slash groups seem irrelevant. However, recognizing the large 
amount of variation among sampling points and the similarity of some of the regressions 
in table 3, the equations in figure 4 are recommended for application to fire modeling. 
To obtain the equations for initial depth shown in figure 4, some of the results in 
table 3 had to be adjusted to age 1 year. 


Merchantable top diameter limits and d.b.h. of trees thinned precommercially 
appear to significantly influence the depth-intercept count relationship. We were 
unable to include a wide range of conditions for these factors; thus, application of 
our results should be restricted to conditions similar to ours. Applicable conditions 
include merchantable top diameters of 4 to 6 inches for pine and 5 to 6 inches for other 
species, and precommercial thinning of trees 2 inches and greater in d.b.-h. 


TABLE 3.--Constrained regressions through aggregated data points relating average bulk depth, y (inches) to average number 
of intercepts of 1/4- to 1-inch debris fuel in two-crossed 2-foot sample planes, x 


2 : Age of : : : 
Type of cutting : Dominant : debris : Sample included for data averaging : Regression results: Figure 
: species : (year) : Year Units : Equation "Ss"; reference 
Ground-lead harvest, 6-in top DF, mix 1 1975 Pp ER ils alt! y= 3.16v¥x 0.86 
Ground-lead harvest, 6-in top 
(slashed after harvest removed) DF 2S 1975 1; 85 9 y = 3.22 87 
Ground-lead harvest, 6-in top 
(slashed after harvest removed) DF 1 1975; 1105) 12 y = 3.18Vx -86 
Ground-lead harvest, 6-in top DF 3 1975 3,5 y = 3.45vx .71 
Ground-lead harvest, 3-in top LP Me 1975 255 245, 219 20s y = 3.03¥x 88 
1976 3, 4 
Ground-lead harvest, 6-in top DF, mix NAARS 1975 2955051059 LO LL LZ, 1S, 14 y = 3.27V¥x 195 3 
Ground-lead harvest, 5-in top LP 1 976° 23,5, +24,.150 y = 4,.51v¥x 84 
Precommercial thinning PP 4 1976 LOS Zs, SS a6 y = 3.56v¥x 92 
Precommercial thinning!/ LP 4 1976 I 535. 6)5. Liz y = 5.22v¥x 98 
Precommercial thinning DF Bt4 975; 974, 1557 16,, 17, 18s 19s 22 y = 4.75vx -90 1 
Precommercial thinning LP, PP 1 1976 Ose 27> 26,0295 OL y = 5.58v¥x 94 
High-lead harvest, 6-in top Mix 1 1976 18, 20,21, 26,35,36,37,38,39 y = 4.60V%x 91 2 
High-lead harvest, 6-in top 
(lopped by harvest crew) Mix 1 1976 Tar OR 22,525 y = 3.60Vx .96 


1/ 


— These units are not typical of current practice, as trees of large d.b.h. were thinned, The relationship 
developed was not used for model purposes for this reason, 


BULK DEPTH (INCH) 


Figure 1.--Aggregated data 
for bulk depth as a func- 
tion of intercept count 
(1/4- to 1-tnch particles 
per 4 feet of transect) 
for precommeretal thinning. 
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Figure 2.--Aggregated data 
for bulk depth as a fune- 
tion of intercept count 
(1/4- to 1-tneh particles 
per 4 feet of transect) 
for high-lead skidding. 
Analysts ts shown in 
table 3. 


Pigure 3.--Aggregated data 
for bulk depth as a fune- 
tton of intercept count 
(1/4- to 1-tneh particles 
per 4 feet of transect) 
for ground-lead skidding. 
Analysts ts shown tn 
table 3. 
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Bulk Depth Versus High Intercept Depth 


Measurement of bulk depth is time consuming and subject to personal interpretation. 
However, the simpler to measure "high intercept depth" (Brown 1974) tends to under- 
estimate fuel bed compactness, resulting in a systematic overestimation of the rate of 
fire spread and reaction intensity by Rothermel's (1972) model (Brown 1972; Bevins 1976; 
Hough and Albini 1978). But because it is rapid, reliable, and widely used, we deter- 
mined the regression relationship between high intercept depth and bulk depth to provide 
a formula that can be used to reduce the high intercept measurements for use in the 
Rothermel fire spread model. The analysis was made on pairs of bulk depth and high 
intercept depth measurements taken at each point. 


Examination of a large number of regressions for various combinations of skidding, 
species, and age revealed no substantive differences between the regression coefficients. 
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Also, all the regression lines passed very nearly through the origin. For example, 
treating all debris of age 2 years and less as one set and all debris older than 2 years 
as another set gave the following relationships: 


AGE <2: y = -0.0245 + 0.661X (r2 = 0.78) 
AGE*S2> 0.303 + 0.602X (x2 = 0.74) 


where 
bulk depth (inches) 
x = high intercept depth (inches). 


“< 
i 


Variation among regression coefficients was restricted as shown by the histogram 
of ratios of bulk depth-to-high intercept depth computed for 118 transects (fig. 5). 
Sixty percent of the ratios were between 0.55 and 0.75. Regression coefficients 
pertaining to the major skidding, species, and age groups are shown in table 4. Vari- 
ability of coefficients within slash groups is comparable to variability among slash 
groups. Thus, the narrow range of variability and probable difficulty in establishing 
significant differences among slash groups seems to warrant application of one rela- 
tionship between bulk depth and high intercept depth to all slash. Combining all 
measurements produced the regression: 


= 0.638X 
(s = 0.76) 


which simply states that bulk depth is 64 percent of high intercept depth. This com- 
pares with 52 percent observed by Bevins (1976) in studying Douglas-fir and hemlock 
slash in Washington. 


Interestingly, the finding of 64 percent is consistent with fire spread verifica- 
tion studies (Brown 1972; Bevins 1976; Hough and Albini 1978) in which fuel bed depths 
that were reduced by factors in the range of 0.6 to 0.7 improved fire spread predictions 
by the Rothermel model. 
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Figure 5.--Htstogram of 
regression coefftctents 
for the constrained 
relattonshtp between 
butk depth versus high 
intercept depth deter- 
mined for tndivtdual 
transects. 
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BULK DEPTH / HIGH INTERCEPT DEPTH 


TABLE 4.--Constrained-regresston coefficients (b = ratto of bulk depth to high-intercept depth) 
for skidding, spectes, and age groups 


: svAgeso£ : b-coefficient 
Type of cutting : Dominant : debris : Areas included : Average Standard 
$ species ; (years) : 5 deviation 
Ground-lead harvest DF 0-1 1975-2,10,11,12,13,14,20 0.610 0.030 
Ground-lead harvest DF 3-4 1975=3,5,,7,8,9 .599 031 
Ground-lead harvest LP 3 1975-23,24,25,26 596 .079 
Ground-lead harvest LP 1 1976-3,4,23,24,30 698 .036 
Precommercial thinning PP 4 1976-10,12,13,15,16 2092 L07 
Precommercial thinning LP 3-4 V9O76=1)5.2'5:5:,'6:, 7 lor. 653 . 066 
Precommercial thinning DF 3-4 1975-4,15,16,17,18, 195,215.22 55175 .092 
Precommercial thinning PP 1 1976-9,11,27,28 sHs .059 
Precommercial thinning LP 1 1976-29,31 .689 . 086 
Skyline harvest Mixed 1 1976-14,18,19,20,21,22,25,26 1673 . 036 
Helicopter harvest Mixed 1 1976-35,36,37,38 -631 .048 
Effect of Lopping 


The effect of lopping was evaluated by combining bulk depths determined before and 
after lopping, as a data pair for each study area and computing a linear regression 
constrained through the origin. The resulting equation is plotted over the scatter 
diagram shown in figure 6. Similarly, the data from all precommercial thinning areas 
were combined for regression (as shown in fig. 7). The quality of the relationships 
expressed by "s'' values of 0.91 and 0.69, respectively, was good. Note that the ratio 
of the regression coefficients in table 3 for lopped and unlopped high-lead harvest 
debris (3.60/4.60 = 0.78) is quite close to our result of 0.83 (fig. 6). 


HAZARD Model Application 


To illustrate how this study has been utilized, sample printouts of fuel inputs 
(fig. 8) and fire behavior predictions (fig. 9) are shown for the slash HAZARD model 
now in use by USDA Forest Service Northern Region. This model affords managers the 
opportunity to assess the fire implications of tree cutting activities before debris is 
put on the ground. The fuel loadings required to run the HAZARD model are generated 
from tree inventories processed through a debris prediction model.° The debris pre- 
diction model provides total potential debris for all trees on a site. Managers then 
tailor the total potential debris to specific cutting prescriptions and submit the 
data for processing by the HAZARD model. To ease checking of the transferred data, 
the HAZARD model prints out the input data (fig. 8). 


5 See footnote 4. 
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Figure 6.--Bulk depth of 
lopped slash as a func- 
tton of bulk depth of 
unlopped slash for all 
timber harvest conditions. 


Figure 7.--Bulk depth of 
lopped slash as a function 
of bulk depth of unlopped 
slash for precommerctal 
thinning. 
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EXAMPLE IN®UT/OUTPFUT FOR HAZARD PROGRAM 

ORIGINAL CONDITION 
POTENTIAL, FIRE BEAAVIOR FOR SEVERELY=-DRY CONDITIONS-= 
(ONE-HOUR TIMELAG FUEL MOISTURE CONTENT = 5S PERCENT) 


QUANTITY 20-F Te SLASH AGE IN NUMBER OF WINTERS 


ASSESSED WIND: 1 3 5 
MI7Hd sess £“seses § (see=— 
@eeoeesnee8 ee @ 0 4 3 l 
HEAD FIRE 5 7 5 3 
SPREAD RATE 10 re 8 4 
CHAINS/HOUR 15 19 13 6 
e@eee2eenea2eoede00@ 20 27 18 8 
@eeoeaeeneetdeeo@ 0 15 ll 6 
PERIMETER 5 24 17 10 
GROWTH RATE 10 36 24 13 
CHAINS/HOUR 15 47 31 16 
e@@eeeeeeaoeaoeeeoe 20 57 37 17 
@eeeeteeeoee08@ 0 2 1 0 
BURNED AREA 5 5 2 1 
AFTER 1 HAOJR 10 11 5 1 
BURN® ACRES 15 20 9 2 
@@eeoesteeoeoeee @ 20 29 12 3 
@eoeoet®oeoeeee @& 0 6 5 3 
HEAD FIRE 5 2 7 5 
FLAME. LENGTH 10 11 8 6 
IN FEET 15 14 10 i 
coccccccccce 20 16 re 8 
e@oeeeteoneatded @ 0 4&7 31 19 
CROWN SCORCH i) 65 4) ee 
HEIGHT IN FT 10 16 43 18 
BY HEAD FIRE 15 82 42 16 
@eeeeteeneeteoe@ 20 85 42 14 
@eeeeeeneeee @ 0 324 176 Bl 
INTENSITY AT 5 600 335 158 
HEAD OF FIRE 10 1081 584 254 
(BTU/SEC/FT) 15 1688 885 358 
e@eeeteooaoeeoed 20 2395 1224 468 


RESIOUAL FIRE FROW BURNOUT OF LARGER FUELS (AND DUFF) 


HEAT LOAD (BTU/SQFT) 16439 15202 14344 
BURN TIME (MINUTES) 69 69 69 


GRDe LEAD HARVEST UTLe STD- 620 IN RUN 14 NOV 77 


Figure 9.--Sample printout of ftre behavior predictions by HAZARD program, based on 
fuel tnputs shown tn figure 8. 
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Fire behavior predictions are furnished in the same format for unlopped and lopped 
debris (as shown in fig. 9). A further discussion of model output and a guide to inter- 
pretation of the fire behavior numbers are in the Northern Region's Users' Guide.® 


Fuel Bed Depth Prediction 


To predict fuel bed depth, predictions of slash in tons per acre are converted to 
number of 1/4- to l-inch intercepts through two manipulations. First, weight per unit 
area is converted to number of 1/4- to l-inch intercepts by species, using a constant 
multiplier that is inversely proportional to the product of wood density of the 1/4- 
to l-inch pieces and their mean square diameter (Brown 1974). Table 5 gives the 
conversion factors for 11 western conifer species. 


Table 5.--Factors for converting from 1/4- to 1-inch size class 
fuel loading to "expected number of tntercepts per 
foot of random transect perpendteular to the ground" 
for 1/4- to 1-tnch stze class limb wood of dtfferent 


spectes 
Species of tree : Factor for tons/acre.—: Factor for lb/ft 
Ponderosa pine 0.4309 91939 
Lodgepole pine 55871 12279 
Western larch 4382 9...55 
Douglas-fir .5608 V2; 22 
Grand fir M295 15.90 
Subalpine fir .5784 12.60 
Western redcedar .6400 13.94 
Western white pine aD298 11.46 
Western hemlock -5950 12.96 
Engelmann spruce . 7089 15.45 
Whitebark pine “258 11.46 


Next, number of intercepts are combined for different species. Equations have been 
given that relate bulk depth to the count of 1/4- to 1l-inch fuel intercepts in two 
crossed, 2-foot vertical planes (this relationship is the same as that for a randomly 
placed single 4-foot vertical plane). The equations for the pines and for other 
conifers are different, so in some cases a method of combining them is necessary. 


The depth prediction equation can be written for slash of one species type as: 


eS ese 
fe) 
where 
65 = initial bulk depth 
xX = expected number of intercepts of 1/4- to l-inch fuel particles in a randomly 


placed 1-foot vertical plane. 


®See footnote 3. 


ey 


The quantity x is proportional to the average loading on the site of 1/4- to l-inch fuel 
pieces. When only pines are present in the fuel bed, or when none are present, this 
formula provides a prediction of fuel bed depth for fire modeling. But when both types 
are present, the mixed species fuel bed depth will contain contributions for both types. 
The model for mixed type fuel bed depths used in the HAZARD model is as follows: 


Let xy expected 1/4- to l-intercept count per foot for pine types, 


and X9 expected 1/4- to l-intercept count per foot for other types. 


If the two types are randomly distributed over the site, then the fraction of the total 
1/4- to l-inch size class loading that is pine type should be (ignoring particle density 
differences) f£), where 


£, = x1/(X) + Xg) 


and the fraction for other types is f,, where 
fo = Xo/ (X, Fao )k v= all eeeeteas 


Now, if the two types were segregated on the site, so that a fraction f; of the 
site area would be covered by pine species and f2 by others, the expected intercept 
count tm the pine-covered area would be x)/f, and the bulk depth in the pine-covered 
area would be 


oe =.2) ¥X,/f£y- 


Similarly, the bulk depth of the other fraction of the area would be 


oe) ES an/x,/f,. 


On such a segregated site, the average bulk depth would be 85° where 


ee (1) (2) 
84 = £155 ot £26 5 : 


We take this latter expression to be the average fuel bed depth for the mixed type sit- 
uation. This equation can be rewritten by substituting for the f and oF values to give 


ie = (a,X} aA AXo)/VXy aig Xo. 


Other factors that affect initial fuel bed depth are reflected in this formula by chang- 
ing the values of a, and a, in accordance with the results discussed earlier (fig. 4). 


Foliage Loss and Settling 


The early effects of aging on debris fuels are settling of the fuel bed and loss of 
foliage and fine twigs to the forest floor. For modeling fire behavior in slash, it is 
necessary to quantify rate of settling and the loss of foliage and fine twigs. To 
operate the HAZARD model, data from Olson and Fahnestock (1955), Fahnestock and Dieterich 
(1962), Steele (1960), Wagener and Offord (1972), and Brown (1970a) were integrated to 
describe loss of material (fig. 10). Not all foliage that drops from the branches was 
excluded from the fuel complex because some foliage remaining in the litter layer and 
suspended as mats in the slash is still available as fuel for a surface fire. 


Settling of slash was modeled as a reduction in depth using table 3, Fahnestock 


and Dieterich (1962), and Kiil (1968). This study, together with others cited, permitted 
construction of a settling model adequate for hazard appraisal (fig. 11). 
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Figure 10.--Reductton of 
quantity of foltage and 
small twigs for modeling 
fire behavior in slash. 
Not all foltage detached 
from limbs ts lost from 
fuel avatlable to spread 
fire. 


Figure 11.--Settling of 


slash expressed as a 
funettion of the ortginal 
depth over ttme for har- 
vesting and precommerctal 
thinntng slash. 
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SUMMARY 


To aid in modeling fire behavior for appraisal of slash hazard, equations were 
developed for predicting fuel bed depth from the loading of 1/4- to 1l-inch fuels. The 
depth equations differed among precommercial thinning, ground-lead harvesting, and 
high-lead harvesting systems and between pines and other conifers. Settling of slash 
with age was weakly discerned; however, our findings, together with others, permitted 
inference of a rough aging model that distinguishes between harvesting and precommercial 
thinning. 


Depth of lopped slash was strongly related to depth of unlopped slash. A strong 
relationship was also developed between the ‘bulk depth," useful in fire behavior 
modeling, and the "high intercept'' depth that is easily measured in the field. The 
influence of merchantable tip diameter on slash depth was evident, although it was not 
quantified. Accuracy of the depth relationships is probably adequate for most tree cut- 
ting activities in western mountainous areas, except perhaps west side slopes along the 
Pacific Coast. Here the large trees and methods of logging may result in depth/load 
relationships different from ours. Depth predictions should be verified in the Pacific 
Coast mountains before modeling fire behavior in slash. 


When predictions of fuel depth are coupled with predictions of slash loading (by 
size class) from tree inventories, the depth models presented here allow prediction of 
fire behavior before the slash is created. Hopefully, this information will aid in 
timber sale planning and in management of slash fuels. 


REFERENCES 


Albini, Frank A. 
1976. Computer-based models of wildland fire behavior: a users' manual. USDA For. 
Serv. Intermt. For. and Range Exp. Stn., 68 p. Ogden, Utah. 
Bevans,, Collin. 
1976. An evaluation of the slash fuel model of the 1972 National Fire Danger Rating 
System. (MoS thesis, Univ. Wash., Seattle. 104 p. 
Brown, James K. 
1970a. Vertical distribution of fuel in spruce/fir logging slash. USDA For. Serv. 
Res. Pap. INT-81, 9 p. Intermt. For. and Range Exp. Stn., Ogden, Utah. 
Brown, James K. 
1970b. Ratios of surface area to volume for common fine fuels. For. Sci. 16:101-105. 
Brown, James K. 
1972. Field test of a rate-of-fire-spread model in slash fuels. USDA For. Serv. Res. 
Pap. INT-116, 24 p. Intermt. For. and Range Exp. Stn., Ogden, Utah. 


2a 


Brown, James K. 

1974. Handbook for inventorying downed woody material. USDA For. Serv. Gen. Tech. 

Rep. INT-16, 24 p. Intermt. For. and Range Exp. Stn., Ogden, Utah. 
Brown, James K. 

1978, Weight and density of crowns of Rocky Mountain conifers. USDA For. Serv. Res. 

Pap. INT-197, 56 p. Intermt. For. and Range Exp. Stn., Ogden, Utah. 
Brown, James K., and Peter J. Roussopoulos. 

1974. Eliminating biases in the planar intersect method for estimating volumes of 

small fuels. For. Sci. 20(4):350-356. 
Byram, George M. 

1959. Combustion of forest fuels. Jn Forest fire control and use, p. 90. Kenneth 

P. Davis, ed. McGraw Hill Book Co., New York. 
Fahnestock, George R., and John H. Dieterich. 

1962. Logging slash flammability after 5 years. USDA For. Serv., Intermt. For. and 

Range Exp. Stn. Res. Pap. 70, 15 p. Ogden, Utah 
Frandsen, Wiliiam H. 
1974. Spatial nonuniformity of slash fuels from line transects. USDA For. Serv., 
North. For. Fire Lab. Study Plan FS-INT-2103-09, 7 p. Missoula, Mont. 
Harmathy, T. Z. 
1972. A new look at compartment fires. Fire Technol. 8(3):196-217, and 8(4):326-351. 
Hough, W., and F. A. Albini. 

1978. Predicting fire behavior in palmetto-gallberry fuel complexes. USDA For. Serv. 

Res. Pap. SE-174, 44 p. Southeast. For. Exp. Stn., Asheville, N.C. 
Krad. AvoD: 

1968. Changes in the physical characteristics and moisture content of pine and 
spruce/fir slash during the first 5 years after logging. For. Res. Lab. Intern. 
Rep. A-14, 49 p. Edmonton, Alberta, Can. 

Olson, D. S., and George R. Fahnestock. 

1955. Logging slash--a study of the problem in Inland Empire forests. Univ. Idaho 

For., Wildl., and Range Exp. Stn. Bull. 1, 48 p. Moscow, Idaho. 
Rothermel, Richard C. 

1972. A mathematical model for predicting fire spread in wildland fuels. USDA For. 

Serv. Res. Pap. INT-115, 40 p. Intermt. For. and Range Exp. Stn., Ogden, Utah. 
Roussopoulos, Peter J., and Von J. Johnson. 

1975. Help in making fuel management decisions. USDA For. Serv. Res. Pap. NC-112, 

16 p. North Cent. For. Exp. Stn., St. Paul, Minn. 
Steele, Robert W. 

1960. Fuel hazard rating of slash on state and private lands in Montana. Mont. For. 

and Conserv. Stn. Bull. 14, 18 p. Univ. Mont., Missoula. 
Thomas, P. H. 
1962. The size of flames from natural fires. Jn Ninth Symp. (Int.) on Combustion 
Proc., p. 844-859. Academic Press, Inc., New York. 
Van Wagner, C. E. 
1973. Height of crown scorch in forest fires. Can. J. For. Res. 3(3)2573-378. 
Wagener, Willis W., and Harold R..Offord. 

1972. Logging slash: its breakdown and decay at two forests in northern California. 
USDA For. Serv. Res. Pap. PSW-83, 11 p. Pac. Southwest For. and Range Exp. Stn., 
Berkeley, Calif. 

Williams, F. A. 

1977. Mechanisms of fire spread. In: Sixteenth Symp. (Int.) on Combustion Proc., 

p. 1281-1294. The Combustion Inst., Pittsburgh, Pa. 


U.S. GOVERNMENT PRINTING OFFICE: 1978-0-777-095-34 


22 


Albini, Frank A., and James K. Brown. 
1978. Predicting slash depth for fire modeling. USDA For. Serv. 
Res. Pap. INT-206, 22 p. Intermountain Forest and Range Ex- 
periment Station, Ogden, Utah 84401. 


Development of equations for predicting fuel bed depth (called "bulk 
depth'"' herein) appropriate for modeling fire behavior in slash is de- 
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Application of the models ina slash hazard appraisal computer program 
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